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REPORT  SUMMARY 


This  document  reports  progress  on  ARPA  contract  DAHCO^t-72-C- 
0001  entitled  "ILLIAC  IV  Applications  Research  at  the  Center  for  Advanced 
Computation,  University  of  Illinois  at  Urbana-Champaign."  The  principal 
objective  of  this  program  is  the  development  and  testing  of  numerical 
techniques  and  software  systems  for  use  of  ILLIAC  IV  over  the  ARPA  Net¬ 
work.  This  is  being  accomplished  through  activities  in  the  following 
areas: 

1.  Development  of  numerical  techniques  suitable  for 
parallel  processing  in  the  following  areas: 

a.  Computational  Methods  in  Linear  Algebra. 

b .  Linear  Programming . 

c.  Ordinary  and  Partial  Differential  Equations. 

d.  Time  Series  Analysis. 

e.  Quadratic  Assignment  Code. 

f.  Graphics. 

2-  Development  of  ARPA  Network  facilities  and  associated 
systems  consisting  of: 

a.  Project  to  interface  the  B6700  to  the  ARPA 
Network. 

b.  ARPA  Network  Terminal  System  (ANTS). 

c.  Center  Graphics  Support . 

d.  Network  Graphics  Effort. 

3.  Development  of  a  large  scale  applications  system,  deeding 
with  input-output  economic  modeling,  utilizing  ILLIAC  IV 
algorithms  and  several  ARPA  Network  facilities. 

k.  ILLIAC  IV  Language  Development. 

5.  ILLIAC  IV  Information  Retrieval  and  Statistical  System 
( IRSS) . 

6.  ILLIAC  IV  Image  Processing. 

In  addition,  education  of  segments  of  the  ILLIAC  IV  user  community  was 
accomplished  through  seminars,  classes,  and  the  development  and  dis¬ 
semination  of  tutorial  materials. 
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1.  ALGORITHM  DEVELOPMENT  GROUP 


1.1  Introduction 


The  main  objective  of  this  group  is  the  development  of  numeri¬ 
cal  techniques  that  are  most  suitable  for  parallel  machines ,  namely  the 
ILLIAC  IV.  The  research  program  can  be  divided  into  the  following  areas: 
(a)  Computational  Methods  in  Linear  Algebra 
(b }  Linear  Programming 

(c)  Ordinary  and  Partial  Differential  Equations 

(d)  Time  Series  Analysis 

(e)  Quadratic  Assignment  Code 

(f)  Graphics 

In  developing  any  parallel  algorithm,  we  first  implement  the  existing 
serial  algorithm  on  the  B6700  (in  Algol).  Once  the  programmer  is  familiar 
with  the  computational  characteristics  of  that  routine,  the  parallel  algo¬ 
rithm  is  then  designed.  In  many  instances,  the  parallal  algorithm  is  sub¬ 
stantially  different  from  the  serial  one,  in  which  case  the  parallel  rou¬ 
tine  is  written  on  the  B67OO  in  Algol  to  check  its  validity,  accuracy, 
rate  of  convergence,  etc.  The  last  step  is  then  writing  Glypnir  or  ASK 
code,  and  finally  debugging  it  on  the  B6700  simulator  (SSK).  Due  to  the 
slow  speed  of  the  simulator,  many  codes  cannot  be  debugged  completely. 

Some  time  consuming  codes  have  been  coded  in  Glyplit  (Glypnir  to  PL/1 
translator  developed  at  RAND)  for  fast  checking  of  their  logic  only. 

The  group  has  adopted  a  verification  process  of  the  algorithm 
in  which  each  member  tries  to  use  and  test  algorithms  developed  by  other 
members.  This  leads  to  the  discovery  of  well-hidden  bugs  and  helps  in 
pointing  out  areas  of  needed  standardization  and  documentation,  as  well 
as  improving  the  codes  being  reviewed.  A  listing  of  the  Center’s  Sub¬ 
routine  library  is  now  available,  {l],  which  identifies  available  de¬ 
bugged  algorithms,  describes  their  function,  the  language  in  which  they 
are  written,  and  the  author. 

1.2  Computational  Methods  in  Linear  Algebra 
1.2.1  Solution  of  Systems  of  Linear  Equations 

(a)  An  ASK  program  was  developed  to  solve  systems  of  linear 
equations  Ax  =  b  for  many  right  hand  sides,  where  A  is  a  banded  positive 
definite  symmetric  matrix.  Work  has  been  completed  on  the  algorithm  to 
handle  matrices  of  band-width  <  64,  and  the  algorithm  for  bandwidth  <  127 
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is  currently  in  the  rabugging  state. 

Enroute  to  the  ASK  algorithm,  two  routines  in  Algol  were  extrac¬ 
ted  from  uilkinson  and  Reinsch’s  Linear  Algebra,  [2],  for  incorporation 
the  group’s  subroutine  library: 

(i)  a  routine  which  decomposes  the  banded  matrix,  and 

(ii)  a  routine  which  given  the  decomposition,  proceeds 
to  determine  the  solution  vectors  for  many  right 
hand  sides . 

Documentation  on  the  Algol  program  is  completed,  and  the  ASK 
programs  will  be  documented  upon  their  completion. 

(b)  A  program  for  the  decomposition  of  a  matrix  into  the  pro- 
?■  of  low«.r  and  upper  triangular  matrices  for  the  solution  of  linear 
equations  or  matrix  inversion,  [2],  has  been  completed  in  ASK  for  non¬ 
core  contained  matrices  (order  -  1000).  The  program  has  been  debugged  as 
far  as  possible  but  the  slow  speed  of  the  simulator  and  lack  of  satis¬ 
factory  input/output  routines  have  precluded  further  checking  until  the 
machine  becomes  available.  A  document  is  forthcoming.  A  Glypnir  pro¬ 
gram  for  solving  complex  systems  of  linear  equations,  order  6U,  by  the 
LU  decomposition  has  also  been  completed,  [2, 3, a]. 

(c)  An  ASK  program  of  the  Householder  triangularization,  [3], 
of  a  matrix  and  back-substition  for  solving  core-contained  .crs  terns  of 
linear  equations  (order  =  300),  is  now  in  the  final  stage  of  debugging 
on  the  SSK  simulator.  Be  routine  consists  of  seven  external  subrou¬ 
tines  and  a  master  program,  all  of  which  are  checked  out  independently. 

The  task  presently  in  progress  involves  combining  individual  subroutines, 
several  at  a  time,  for  checking  on  SSK  without  exceeding  its  capabilities. 

Householder’s  Modification  method  for  matrix  inversion  has 
also  been  coded  and  debugged  in  ASK  but  only  for  matrices  of  order  6h  or 
less  since  the  number  of  operations  is  higher  than  that  of  the  LU  decoxqpo- 
sition  or  Householder  triangularization.  The  method  seems  to  be  suit¬ 
able  for  sparse  matrices.  A  document  describing  both  Householder's  tri¬ 
angularization  and  the  Modification  methods  is  forthcoming. 

(d)  The  conjugate  gradient  method  for  solving  linear  equa¬ 
tions  Ax  =  b,  [2],  has  been  coded  in  Algol  to  gain  experience  in  the 
technique.  Glypnir  coding  is  proceeding  but  has  been  somewhat  hampered 
by  lack  of  facilities  in  the  language,  namely  the  inability  to  call  the 
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subroutines  as  parameters.  This  algorithm  requires  as  an  input  a  rou¬ 
tine  to  compute  Ay  and  A^y  for  a  given  y  rather  than  an  explicit  repre¬ 
sentation  of  the  matrix  A.  This  should  prove  to  be  very  fast  in  cases 
where  A  is  sparse  and  of  regular  structure. 

1.2,2  The  Algebraic  Eigenvalue  Problem 

(a)  Two  ASK  programs  for  -narallel  Jacobi  and  Jacobi-like  al¬ 
gorithms,  for  real  symmetric  and  real  nonsymnetric  matrices  of  order  6h 
or  less,  have  been  completed,  and  a  document  has  been  published  [5].  An 
improved  parallel  Jacobi  algorithm  has  been  developed  [6],  and  an  ASK 
program  hoc  been  completed  and  documented  [y].  Several  experimental 
Algol  and  Glyplit  prograrjs  have  been  written  to  explore  the  possibility 
of  improving  the  convergence  of  already  published  Jacobi-like  algorithm 
[5,6]. 

(b)  The  QR  algorithm  for  finding  the  eigenvalues  of  symmetric 
tridiagonal  matrices  has  been  coded  in  Glypnir. 

The  QL  algorithm  [2]  for  finding  the  eigenvalues  and  eigen¬ 
vectors  of  symmetric  tridiagonal  matrices  has  been  coded  and  tested  in 
Algol  and  being  programmed  in  Glypnir.  A  generalization  of  the  bisection 
method  together  with  inverse  iteration  have  been  programmed  and  tested  in 
Algol.  Both  routines  are  being  programmed  in  Glypnir  to  produce  a  paral¬ 
lel  algorithm, [8],  with  higher  efficiency  than  the  inherently  serial  QL 
algorithm. 

(c)  The  conjugate  gradient  algorithm  for  finding  the  largest 
and  smallest  eigenvalues  and  the  corresponding  eigenvectors  of  the 
generalized  eigenvalue  problem  Ax  =  Bx  where  A  and  B  are  symmetric  posi¬ 
tive  definite  matrices,  [93,  has  been  coded  and  tested  in  Algol  and  will 
be  coded  in  Glypnir  with  high  efficiency  of  parallelism.  As  mentioned 
earlier,  Az  and  Bz  can  be  generated  implicitly  for  any  vector  z,  when  A 
and  B  structured,  sparse  matrices  (which  is  the  case  in  most  practical 
applications) . 

(d)  The  accelerated  power  method  for  finding  the  largest 
eigenvalue  and  the  corresponding  eigenvector  has  been  coded  and  tested 
in  Algol  and  Glyplit.  This  method,  mentioned  by  Wilkinson  [33  (ad  hoc 
shift  of  origin)  for  roan-machine  interaction,  was  modified  so  that  it 
would  not  require  interaction  unsuited  for  the  ILLJAC  IV.  A  variation 
of  the  above  power  method  now  being  tested,  where  the  norm  is  used 
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instead  of  the  norm,  is  just  as  fast  on  the  ILLIAC  IV  and  more  pre¬ 
dictable,  hence  more  suited  to  unmanned  operations. 

(e)  The  QR  algorithm,  with  origin  shift  for  the  solution  of 
eigenvalue  problems,  [2],  has  been  coded  in  ASK  for  non  core-contained 
matrices  and  is  being  debugged.  The  reduction  to  upper  Hessenberg  form 
using  Householder's  transformations,  [2],  of  non-core- contained  matrices 
is  almost  totally  debugged. 

An  ASK  program  for  reducing  core-contained  real  matrices  to 
the  upper-Hessenberg  form  using  elementary  stabilized  transformation  [2] 
is  completed.  Two  documents  containing  ASK  codes  for  the  QR  algorithm 
and  Householder  reduction  to  upper-Hessenberg  form  for  matrices  of  order 
64  or  less  are  available  [10,11]. 

(f)  An  algorithm  for  establishing  error  bounds  on  the  com¬ 
plete  eigensystem  of  core-contained  (order  -  90)  diagonable  matrices  [12] 
has  been  coded  in  ASK  and  is  in  tht  debugging  stage.  This  should  reveal 
some  of  the  favorable  aspects  of  parallel  computation. 

(g)  A  numerical  method  for  obtaining  the  eigenvector  solution 
of  the  non-linear  matrix  Riccati  equation  has  been  programmed  and  fully 
tested  in  Glypnir,  and  a  complete  document  is  available  [3j. 

(h)  Work,  has  begun  only  recently  on  the  development  of  an 
eigenvalue  algorithm  for  obtaining  all  eigenvalues  and  vectors  of  sparse 
matrices .  As  a  first  step  such  an  algorithm  was  developed  for  matrices 
whose  Gerschgorin  disks  are  all  isolat"*.  The  resulting  parallel  algo¬ 
rithm  has  for  its  basis  the  work  of  Varga  and  his  associates  on  the  best 
isolated  Gerschgorin  disks  [13,16] .  The  Algol  program  has  been  exten¬ 
sively  tested;  results  obtained  compare  favorably  in  accuracy  with  the 
QR  algorithm.  Work  continues  on  the  general  case  where  the  Gerschgorin 
disks  cannot  be  isolated  by  diagonal  similarity  transformations. 

1.3  Linear  Programming 

(a)  The  goal  of  the  initial  LP  codes  is  a  solution  capacity 
of  l6,0G0  row  problems  with  a  potential  for  expansion  to  64,000  row  prob¬ 
lems.  Codes  initially  available  will  be  PRIMAL,  a  modification  of  the 
revised  simplex  method  utilizing  a  product  form  of  the  inverse,  and 
INVERT,  a  reinversion  program  based  on  the  Helleraan-Rarick  preassigned 
pivot  procedures.  Supporting  these  computational  portions  are  the  data 
massaging  SETUP  and  SETD0WH  codes  which  create  initial  files  and  produce 


result-files,  respectively.  The  LOADER,  a  B6500  program,  accepts  the  data  \ 

in  MPS/36O  form  and  after  sorting  passes  this  data  to  the  ILLIAC  IV  disk. 

SETUP  then  utilizes  the  ILLIAC  to  convert  the  values  to  suitable  word  \ 

•  i 

formats,  classify  columns  and  rows,  scale  these  columns  and  rows,  and 
emit  those  files  necessary  for  PRIMAL  and  INVERT.  SETDOWH  picks  up  at 
the  termination  of  the  PRIMAL-INVERT  symbiosis  and  rescales  and  edits  the 
output,  still  using  ILLIAC.  An  UNLOADER  will  eventually  be  written  to 
replace  the  internal  (numeric)  with  the  original  external  (alphanumeric) 
labels  and  to  prepare  printer  output. 

'Hie  Linear  Programming  group,  presently  staffed  by  fewer  then 
two  full-time  staff  members,  progressed  steadily  through  INVERT  coding. 

INVERT  has  been  divided  into  three  sections: 

(i)  the  creation  of  an  implicit  in-core  matrix 
representing  the  basis  vectors, 

(ii)  the  identification  of  the  pivot  sequence  for 
the  new  Product  Form  of  the  Inverse  (PFl), 

(iii)  the  creation  of  the  ETA  files  representing 
the  PFI. 

The  first  two  sections  have  been  coded  and  simulated;  while 
their  lack  of  parallelism  makes  them  less  than  satisfactory,  they  are  the 
best  algorithms  identified  so  far.  The  third  section  appears  to  make  very 
good  use  of  the  parallelism  of  the  machine  and  should  be  simulated  by 
the  end  of  the  present  quarter. 

PRIMAL  has  been  coded  and  debugged  as  fas  as  the  simulator 
could  allow.  The  major  revision  being  considered  is  the  inclusion  of 
debugging  and  diagnostic  routines .  A  lack  of  response  to  our  letters  re¬ 
questing  information  from  NASA/AMES  Institute  for  Advanced  Computation 
(ascribed  to  a  shortage  of  personnel  at  IAC)  has  postponed  a  rewrite  of 
our  I/O  utilities  and  the  development  of  these  diagnostic  routines.  In¬ 
formation  obtained  at  the  March  1972  user's  conference  at  Monterey  suggests 
that  ILLIAC  IV  will  not  begin  to  support  LP's  I/O  requirements  for  some 
time. 

The  SETUP  routines  are  nearly  completed  and,  with  the  current 
resolution  of  the  last  section  of  INVERT'S  flowcharts,  should  be  quickly 
disposed  of.  Documentation  of  all  coded  routines  are  available  and  will 
be  combined  into  one  document  describing  the  whole  system  once  unresolved 
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issues  about  ILLIAC  IV  are  settled. 


lb)  Preliminary  investigation  of  methods  to  reduce  storage 

and  computation  time  for  large  Linear  Programming  problems  have  begun. 

This  would  be  achieved  by  performing  all  matrix  operations  (e.g. ,  cal- 
T 

culating  Ax  and  A  y  given  x,y)  implicitly  without  ever  generating  and 
storing  the  matrix  A.  Used  with  a  matrix  generator,  significant  savings 
may  be  obtained. 

1.1»  Ordinary  and  Partial  Differential  Equations 

(a)  A  preliminary  study  for  an  ordinary  differential  equa¬ 
tions  solver  for  ILLIAC  IV  has  been  completed  and  a  document  is  avail¬ 
able  [17]. 

(b)  Two  Glypnir  codes  have  been  completed  for  the  alternat¬ 
ing  direction  implicit  iterative  technique  (ADI,  6U  rows),  and  the  block- 
Jacobi  method  (128  rows)  for  the  solution  of  second-order  elliptic  par¬ 
tial  differential  equations.  A  document  describing  the  block-Jacobi 
method  on  the  ILLIAC  IV  is  forthcoming. 

( c)  A  Glyplit  program  has  been  written  which  performs  simul¬ 


taneous  fast  Fourier  transforms  of  real  functions  in  all  enabled  PE's. 

At  present,  the  program  is  ccre-contained.  Elis  program  may  be  useful 
when  solving  partial  differential  equations  by  means  of  a  change  of 
variables  that  can  be  carried  out  by  the  fast  Fourier  transform  [l8] . 

(d)  A  Glypnir  program  implementing  J.  Son's  (Department  of 
Chemical  Engineering,  University  of  Illinois  at  Urb ana-Champaign)  solu¬ 
tion  to  the  problem  of  flow  around  a  cylinder  has  been  coded  and  debugged 
on  the  simulator.  It  Is  estimated  that  a  50%  reduction  in  run  time  can 
be  achieved  by  coding  in  ASK  the  section  which  solves  Poisson's  equation 
iteratively. 

Direct  methods  for  solving  Poisson's  equation  appear  to  offer 
additional  reductions  in  run  time.  Hockney’s  direct  method  has  been 
implemented  in  serial  form  in  Algol.  This  algorithm  vas  compared  to  a 
serial  SOR  and  found  to  be  about  twice  as  fast  for  those  cores  tested. 

A  direct  algorithm  for  solving  Poisson's  equation  in  an 
irregular  domain  using  a  Hockney  type  formulation  is  being  developed.  A 
serial  iterative  algorithm  has  already  been  developed  in  Algol  for  these 
irregular  domains. 

Various  boundary  conditions  are  also  being  tested  in  an 


attempt  to  tetter  model  the  physics  of  this  fluid  flow  problem. 

(e)  A  computational  technique  for  high  subsonic  compressible 
inviscid  flow  past  a  circular  cylinder  has  been  implemented  serially  on 
the  CDC  l6o4  computer;  a  Glyplit  code  will  be  written  soon.  Boundary 
conditions  are  chosen  for  maximum  utilization  of  the  processing  elements. 
Procedures  for  varying  grid  sizes  without  changing  the  PE  instruction 
stream  have  been  developed.  Ime  calculated  results  compare  well  with 
those  obtained  by  Holt  and  Kasson  [l8]  by  the  method  of  integral  rela¬ 
tions.  A  document  has  also  been  completed  [19]  in  which  the  relevant 
and  practical  criteria,  required  for  the  numerical  calculation  cf  differ¬ 
ent  types  of  gas  dynamic  flow  fields  are  discussed. 

(f)  Glypnir  programs  have  been  written  to  solve  a  set  of  differ 
ential  equations  used  in  atmospheric  dynamics  to  describe  fluid  convec¬ 
tive  motion  induced  by  the  buoyancy  force  and  heat  transfer  associated 
with  the  motion.  Time  comparisons  have  been  made  between  these  programs 
and  CDC  6 4 00  FORTRAN  programs  for  a  64  x  17  mesh.  The  ratio  is  about 
1:110  when  successive  overrelaxation  method  is  used  in  solving  Poisson’s 
equation.  To  improve  this  ratio,  the  kernel  of  this  program  was  pro¬ 
grammed  in  ASK.  The  ratio  is  thereby  improved  to  1:220  as  the  result  of 
storing  constants  in  ADB's,  making  more  efficient  address  calculations, 
and  increasing  the  amount  of  overlap  between  PE  and  PEM.  Detailed  analy¬ 
sis  of  the  kernel  of  Glypnir  and  partially  ASK-coded  programs  have  been 
performed. 

In  solving  Poisson’s  equation,  the  efficiencies  of  three 
popular  iterative  methods  (SOR,  SLOR  and  ADI)  have  been  compared  on 
ILLIAC  IV.  The  storage  schemes  and  methods  of  implementation  which 
maximize  the  efficiency  of  these  algorithms  on  ILLIAC  IV  have  been  in¬ 
vestigated.  The  implementation  of  the  direct  method  in  solving  Poisson's 
equations  on  rectangular  regions  with  various  boundary  conditions  is 
under  vTay.  Detailed  reports  of  this  work  will  be  available  next  quarter. 
These  studies  on  the  efficiency  of  ILLIAC  IV  in  atmosphere c  dynamics 
calculations  axe  jointly  supported  by  ARPA  and  the  National  Science 
Foundation. 

(g)  Algorithms  for  the  "identification"  of  non-linear  differ¬ 
ential  equations  have  been  developed.  Such  algorithms  will  allow  deter¬ 
mination  of  certain  sets  of  differential  equations  that  satisfy  a  set  of 


data  provided  as  input.  These  routines  have  teen  implemented  on  the 
B67OO  in  Algol;  however,  many  features  are  suitable  for  parallel 
machines.  A  document  will  soon  be  published.  We  believe  that  the 
developed  package  -will  be  of  interest  to  many  of  the  ARPA  community  who 
are  involved  in  mathematical  model  building.  If  the  response  is  favor¬ 
able,  work  will  proceed  on  implementing  the  package  on  ILLIAC  IV.  This 
is  the  first  of  three  studies  begun  at  the  Center  for  mathematical 
model  building.  The  other  two  studies  are  stochastic  model  building  and 
identification  of  linear  and  periodic  dynamic  systems  which  are  described 
in  Section  1.5. 

1*5  Tine-Series  Analysis 

(a)  A  Fast  Fourier  Transform  (FFT)  subroutine  has  been  writ¬ 
ten  in  ASK,  and  fully  debugged  and  tested.  Timing  simulation  results 
have  shown  rhat  this  program  will  perform  Fourier  transforms  on  ILLIAC 
TV  with  good  efficiency.  A  descriptive  document  and  a  user’s  manual  are 
available  [21]. 

(b)  Glypnir  codes  for  a  one-dimensional  FFT,  autocorrelation, 
and  computation  of  the  impulse  response  of  a  hand-pass  filter  are  fully 
debugged  and  tested.  A  document  is  forthcoming. 

(c)  A  document  has  been  completed  [22]  describing  a  B6700 
package  developed  to  determine  the  stochastic  model  used  to  describe 
observed  data  (which  is  in  the  form  of  time  series),  determine  the  order 
of  the  model,  and  estimate  the  parameters  using  a  maximum  likelihood. 

The  feasibility  of  using  parallel  computations  in  this  study,  especially 
for  the  estimation  algorithm,  is  being  considered. 

(d)  Algorithms  have  been  developed  for  the  identification  of 
linear  time-invariant  and  periodic  (period  is  not  known  a  priori)  dyna¬ 
mic  systems  using  only  the  observed  data  (no  control  inputs).  A  B6700 
Algol  package  is  now  being  completed.  The  obtained  results  are  very 
satisfactory.  Many  parts  of  the  developed  algorithms  are  suitable  for  a 
parallel  computer,  especially  since  the  identification  of  a  periodic 
system  is  reduced  to  the  identification  of  many  (64)  linear  system. 

Again,  we  believe  this  study  to  be  valuable  for  researchers  involved  in 
mathematical  model  building  (Economics,  Biology,  Engineering,  etc.). 

1.6  QUASCO:  An  ILLIAC  IV  Quadratic  Assignment  Code 

As  part  of  the  ILLIAC  IV  algorithm  development  effort  (and 
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in  support  of  CAC  hardware  design  research)  a  modest  effort  is  underway 
to  formulate  and  code  efficient  parallel  procedures  for  treating  the 
quadratic  assignment  problem.  The  results  of  these  investigations  to 
date,  though  incomplete ,  seem  encouraging. 

The  quadratic  assignment  problem  requires  a  one-to-one  match¬ 
ing  of  the  elements  of  two  sets  such  that  a  quadratic  objective  function 
is  optimized.  For  example,  in  the  area  of  computer  hardware  design  it  is 
often  desirable  to  assign  discreet  positions  on  a  wirins  board  to  some 
number  of  interconnected  components  in  such  a  manner  that  the  total 
length  of  connections  between  components  is  minimal.  Other  applications 
of  the  model  relate  to  optimal  industrial  plant  layout,  hospital  design, 
office  space  allocation,  urban  planning,  etc. 

A  possible  notation  for  the  problem  (in  simplified  form)  is: 


min  I  I 

P 

i=l  j=l 


“p(i),  p(j) 


4i.J 


where: 


[Cp(l)  o^j)  =  annxn  matrix  expressing  the  cost  per  urdt 
distance  between  any  two  elements  of  the  first  set  ( components ,  plant 
functions,  departments ,  offices,  activities,  etc.). 


[D.  ,]  =  an  n  x  n  matrix  expressing  the  effective  distance 

x,j 


between  any  two  elements  of  the  second  set  (positions  on  a  printed  cir¬ 
cuit  board,  spaces  within  a  plant,  hospital  or  office  tower,  city 
blocks,  etc.). 

p  =  a  permutation  of  indices  of  the  elements  of  the  first  set 
corresponding  to  a  one-to-one  matching  to  the  elements  of  the  second  set. 

While  integer  programming  formulations  of  the  problem  are 
possible,  the  effect  of  the  resulting  quadratic  objective  function  and 
integer  constraints  is  always  to  create  non-convexities  in  the  solution 
space  and  hence  numerous  locally-optimal  solutions .  It  may  indeed  be 
impossible  to  compute  global  optima  for  large  quadratic  assignment  prob¬ 
lems,  even  on  ILLIAC  IV. 

Thus,  our  attack  has  focused  more  on  the  development  of  paral¬ 
lel  procedures  for  approximating  optimal  solutions.  Since  these  methods 
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employ  hlll-clinibing  programming  strategies,  the  solutions  depend  to 
some  extend  on  starting  points  (initial  feasible  assignments  derived  by 
chance  or  design) .  Given  such  algorithms  and  an  appropriate  sampling 
strategy,  however,  the  speed  at  which  ILLIAC  IV  will  be  able  to  generate 
and  rank  alternative  approximations  will  afford  more  confidence  in  the 
calculation  and  employment  of  best  approximations  to  the  optimal  solution. 

An  existing  quadratic  assignment  code  developed  by  a  staff 
member  was  re-written  in  Algol  for  the  E6700  in  order  to  realize  effi¬ 
cient  computer  usage  and  inexpensive  costs.  Extersive  improvements  were 
made  to  the  algorithm;  however,  it  still  remains  in  a  state  of  dynamic 
change.  Test  results  thus  far  indicate  that  QTJASCO  has  achieved  better 
results  than  other  known  algorithms.  Work  will  proceed  on  an  ILLIAC  XV 
routine  when  this  heuristic  algorithm  becomes  better  formulated. 

A  report  is  forthcoming  in  which  are  presented  several  ideas 
on  the  quadratic  assignment  problems  and  the  approximation  of  one  matrix 
by  another  of  lower  rank. 

1-7  Graphics  Algorithm  Development  for  ILLIAC  TV 


1.7-1  Surface  Presentation  of  a  Function  F(X,Y) 

A  basic  ASK  routine  was  completed  for  the  perspective  pre¬ 
sentation  of  a  graphical  image  of  a  surface  represented  by  the  function 
of  two  variables,  X  and  Y.  The  design  and  implementation  of  this  algo¬ 
rithm,  however,  is  based  on  the  still  not  rigorously  defined  input/output 
and  disk  access  aspects  of  the  ILLIAC  IV  operating  system.  Further  work 
in  sending  the  graphics  data  stream  of  this  algorithm  from  ILLIAC  IV  intc 
its  control  system  uas  been  suspended  pending  the  suitable  definition  of 
the  control  j  ystem  by  Ames .  As  soon  as  the  new  I/O  system  to  ILLIAC  IV 
is  defined,  and  the  programming  conventions  established,  this  algorithm 
will  be  remolded  to  fit  the  new  concepts . 

1.7.2  Contour  Plotting  Algorithm  for  ILLIAC  IV 

A  new  contour  plotting  algorithm  for  the  M  x  II  array  case, 
using  real  variables,  was  developed  using  SSK.  The  new  implementation 
is  a  single  program  specifically  designed  to  take  advantage  of  parallel 
processing,  both  in  the  calculation  of  the  contour  data  and  in  the  re¬ 
sulting  production  of  graphical  images. 

The  current  version  of  the  algorithm  assumes  the  case  of  an 
HxH  array  of  data  which  is  entirely  core  contained.  Further  work  on 
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this  algorithm  development  was  halted  pending  complete  definition  of  the 
input/output  system  by  Ames . 

When  work  resumes,  the  case  of  an  M  x  N  array  not  core  con¬ 
tained  (partially  or  wholly  on  disk)  will  be  completed  and  the  input/ 
output  portions  of  the  algorithm  adapted  to  the  new  conventions  for  the 
Ames  system. 

1.8  Miscellaneous 

(a)  Two  algorithms  for  finding  the  roots  of  polynomials  with 
real  coefficients  have  been  programmed  in  Algol-Root-Squaring  and  an 
algorithm  for  finding  the  real  root,s  using  the  Sturm  sequence  property. 

A  parallel  algorithm  is  now  being  designed  to  find  the  real  roots  of 
real  polynomials  using  the  Sturm  sequence  property. 

(b)  A  short  experiment  in  8-fold  precision  on  modular  arith¬ 
metic  was  carried  out  using  eight  prime  numbers  furnished  by  the  DeKalb 
Foundation  for  Humber  Theory. 

(c)  The  block  triangular  decosposition  of  large  matrices  for 
eigenvalue  problems  by  means  of  the  transformation  of  BtPBP"1,  where  P 
is  a  permutation  matrix,  has  been  considered.  The  method  used  finds  the 
transitive  closvre  of  the  digraph  G  associated  with  B. 

An  Algol  program  was  written  and  tested  tc  put  a  b8  x  ma¬ 
trix  in  block  triangular  form,  making  use  of  the  B6T0Q’s  ability  to  lo¬ 
gically  find  the  "and"  of  two  8-bit  words  in  ore  instruction.  Hien  the 
program  was  extended  to  n  x  n  matrices  using  the  same  method.  The  pro¬ 
gram  works,  and  is  now  being  documented. 

The  central  part  of  the  above  program  is  the  calculation  of 
the  transitive  closure  of  a  digraph.  An  ASK  code  to  do  this  for  a  64  x 
6k  matrix  on  ILLIAC  IV  using  Warshall’s  algorithm  [22]  has  been  written; 
it  will  be  extended  to  a  6Un  x  64n  matrix.  Hie  algorithm  is  extremely 
fast  on  ILLIAC  IV  compared  to  conventional  machines. 

(d)  Other  Algol  routines  added  to  the  group's  library  were: 

-  Householder’s  reduction  of  a  symmetric  matrix  to  tri¬ 
diagonal  form. 

-  Hewton's  method  to  obtain  a  root  of  f(x)  =  0. 

-  Wegstein’s  secant  method  to  obtain  a  loot  of  x  =  f(x). 

Hie  above  secant  method  was  programmed  in  ASK  to  run  in  a 

serial  mode  on  ILLIAC  TV. 
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2.  COMPUTER  SYSTEMS  AND  SERVICES  GROUP 

2.1  Introduction 

The  Computer  Systems  and  Services  group  is  responsible  for 
developing  local  basic  computer  systems  and  services  to  serve  the  needs 
of  the  Center  staff  and  user  community.  There  are  four  project  areas  to 
be  reported  on  at  this  time: 

(a)  Project  to  Interface  che  B67OO  to  the 
ARPA  Network 

(b)  ARPA  Network  Terminal  System  (ANTS) 

(c)  Center  Graphics  Support 

(d)  Network  Graphics  Effort 


During  1971*  Burroughs  developed  two  full- capabi li ty  inter¬ 
faces  for  interfacing  a  Burroughs  B6700  computer  system  to  the  ARPA  via 
an  Interface  Message  Processor.  Two  interfaces  were  completely  con¬ 
structed  and  one  of  the  interfaces  was  completely  debugged. 

Software  drivers  for  these  interfaces  were  completed  and  de¬ 
bugged  during  the  Fall  of  1971.  This  included  a  series  of  patches  for 
the  MCP,  as  well  as  the  device  driver  intrinsics  for  the  MCP.  In  addi¬ 
tion,  a  full-size  Network  control  program  was  90?  coded  and  partially 
debugged. 

When  it  was  decided  that  Arnes  would  not  interface  their  B6?00 
to  the  Network,  and  when  it  appeared  that  the  University  of  California 
at  San  Diego  would  write  their  own  network  control  program,  all  effort 
was  stopped  on  the  project.  The  hardware  portion  of  the  project  was 
completed  and  the  NCP  has  been  left  at  the  50?  completion  level. 

It  is  anticipated  that,  in  the  near  future,  UCSD  will  inter¬ 
face  their  system  directly  to  the  Network  and  will  use  one  of  the  hard¬ 
ware  interfaces.  They  will  also  incorporate  our  MCP  patches  and  intrin¬ 
sics  into  their  system.  The  Center  will  aid  UCSD  in  developing  their 
NCP,  but  will  discontinue  an  independent  effort  in  this  area. 

2.3  ARPA  Network  Activities 

2,3.1  ARPA  Network  Terminal  System  (ANTS)  Development 

At  the  beginning  of  the  reporting  period,  development  of  the 
ARPA  Network  Terminal  System  (ANTS)  was  seriously  affected  by  failure  of 
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the  B67OO  performance  to  meet  expectations.  Progress  was  slow,  erratic, 
"but  forward.  Two  full-time  staff  members  have  participated  in  the  de¬ 
velopment  of  ANTS  -  one  devoted  to  developing  the  high-level  language 
compiler  PEESPOL  in  which  the  system  is  written,  while  the  other  was 
responsible  for  the  ANTS  implementation. 

On  September  1,  1971,  ANTS  achieved  Network  operational  sta¬ 
tus  at  the  level  of  terminal  access  to  the  Network.  Since  September  1, 
efforts  have  been  directed  to  consolidating  the  design  and  implementa¬ 
tion  of  the  full  basic  system  in  order  to  facilitate  the  system  func¬ 
tions  described  below. 

As  of  March  1,  the  basic  system  nodules  have  been  completed 
and  the  overall  system  structure  set.  Future  efforts  will  be  directed 
at  implementing  all  applicable  Network  protocol  operations,  such  as  the 
low-level  graphics  protocol,  data  and  file  transfer  protocols,  extended 
terminal  access  protocols,  including  all  required  device  drivers  for 
peripheral  devices,  DECtapes,  magnetic  tapes,  Gould  printer,  card  reader, 
disk  and  graphics  displays .  As  of  April  1,  ANTS  provides  remote  job 
entxy  access  to  the  UCLA  IBM  360/91. 

2.3.2  AEPA  Network  Usage 

With  operational  status  of  ANTS  on  September  1,  1971,  the 
Center  acquired  minimal  capability  for  connecting  to  the  Network.  It 
was  quickly  established,  however,  that  the  Network  was  in  an  early  stage 
and  service  capabilities  were  undeveloped.  Since  September,  experience 
on  the  Network  has  been  primarily  limi'ced  to  attempts  to  use  the  BB&N 
service  site  in  Cambridge  and  the  Network  Information  Center  at  Stanford 
Research  Institute. 

More  Network  sites  have  become  available  during  the  period. 
With  the  recent  addition  of  remote  job  entry  capability  to  UCLA,  Center 
usag^  of  the  Network  will  rise  through  the  use  of  the  GLYPLIT  translator 
on  the  360/91  and  through  the  use  of  the  number-crunching  capabilities 
for  running  large  codes.  Additionally,  the  Center  is  aiding  the  Univer¬ 
sity  of  California  at  San  Biego  in  connecting  their  Burroughs  E6700  to 
the  Network  by  early  June.  This  will  enable  us  to  transfer  work  now 
being  done  on  our  B67OO  to  San  Diego  via  the  Network.  The  Center's  B67OO 
lease  arrangement  terminates  on  June  30,  1972. 
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2.3.3  Further  Installations  of  ANTS  Systems  on  the  Network 

During  the  Fall  of  1971,  the  National  Bureau  of  Standards 
(NBS)  Information  Processing  Techniques  group  acquired  and  installed  an 
ANTS  system.  Due  to  the  hardware  problems,  both  in  the  Network  TIP  and 
the  NBS  PDP-11  system,  this  site  has  not  yet  reached  operational  status. 

It  is  expected  in  the  next  quarter  that  installation  of  the  site  will  be 
complete  as  well  as  installation  of  the  proper  drivers  to  handle  the 
DCT  2000  integrated  remote  job  entiy  terminal  which  IBS  already  has  ac¬ 
quired  and  through  vdiich  they  intend  to  access  the  3o0/91  at  UCLA. 

BBS  intends  to  use  their  version  of  ANTS  to  perform  experiments  and 
gather  data  on  system  performance,  on  the  use  of  terminals  to  various 
ARPA  Network  sites,  and  on  activities  between  systems  on  the  Network. 

They  will  be  modifying  their  operating  version  of  ANTS  through  the  PEESPOL 
compiler  on  the  UCSD  B6700  to  provide  them  with  necessary  input /output 
and  data  collection  capabilities  to  facilitate  these  experiments. 

Additionally,  it' is  expected  that  two  more  installations  of 
the  ANTS  systems  will  he  made  at  Ft.  Belvoir,  Va.  and  Aberdeen,  Kd. ,  both 
sites  of  the  Army  Material  Command  (AMC).  Ft.  Belvoir  will  use  ANTS 
to  interface  remote  job  entry  terminals  to  the  Network  and  to  connect  a 
CDC  66 00  to  the  Network-  Aberdeen  will  obtain  remote  job  entry  access 
to  the  Network  through  ANTS  to  permit  their  AMC  people  to  use  the  Ft. 
Belvoir  CDC  6600  and  other  Network  resources. 

2.^  Graphics  Suooort  for  Center  Projects  on  the  Bo700 


The  development  of  a  graphics  package  for  the  B6700  for  use 
by  the  Center  personnel  was  supported  for  two  reasons:  l)  there  were  a 
number  of  graphics  devices  on  the  PDP-11  which  we  wished  to  m2ke  avail¬ 
able  to  people  losing  the  B6?00;  and  2)  it  was  expected  that  this  exper¬ 
ience  would  aid  in  development  of  graphics  packages  for  ILIIAC  IV  and 
any  future  Network  systems  in  which  we  participated. 

The  connection  of  the  PDP-11  to  the  B6700  was  undertaken  by 
writing  a  special  operating  system  for  the  PDP-11  which  attached  it, 
via  an  lfiOO-baud  telephone  line,  to  the  £6700  and  made  the  PDP-11  peri¬ 
pherals  look  like  terminals  operating  from  the  B6700  Datacoa  processor. 

Further  efforts  were  then  directed  at  providing  a  basic  gra¬ 
phics  package  with  all  rudimentary  graphics  operations  such  as  line  draw¬ 
ing,  point  drawing,  lettering  and  scan  graphics  for  faster  scan 
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conversion  to  the  Gould  electrostatic  printer. 

With  the  completion  of  this  package  in  the  Fall  of  1971, 
efforts  were  then  extended  to  the  development  of  basic  programming  mod¬ 
ules  for  plotting  axes,  scaling  data,  plotting  lines  from  arrays ,  etc. 

In  addition,  a  driving  package  was  developed  which  allows  one  to  inter¬ 
actively  drive  these  programs  to  produce  graphical  imagery. 

The  decision  was  made  to  cease  further  development  of  B67OO 
graphics  in  January  and  to  concentrate  <.  fforts  on  developing  capabilities 
to  the  Network  and  in  the  use  of  other  Network  host  sites  for  Center  per¬ 
sonnel.  Prior  to  this  time,  much  experience  was  gained  on  the  B6700  and 
some  initial  design  on  a  generalized  graphics  data  structure  was  comple¬ 
ted  and  tested. 

2. 5  Network  Graphics  Efforts 

2.5.1  Network  Graphics  Protocol 

Several  staff  members  are  participating  in  the  Network  gra¬ 
phics  group  in  an  attempt  to  establish  protocols  for  the  transmission  of 
graphical  information  across  the  Network.  During  the  first  quarter  of 
1971,  agreement  was  reached  cn  first  level  graphics  protocol  whereby  gra¬ 
phical  images  were  sent  from  the  serving  site  tc  the  using  site  in  a  very 
basic  interpretive  format  allowing  conversion  into  device  specific  for¬ 
mat  to  such  devices  as  the  ARDS .  Tektronix,  Computek  and  IMLAC  graphics 
display  consoles.  Future  efforts  will  expand  this  interpretive  language 
to  provide  for  subroutining  and  dynamic  interpretive  display  of  subpic¬ 
tures. 

2.5.2  Laboratory  for  Atmospheric  Research  (LAR)  Support 

With  the  advent  of  remote  job  entry  to  the  UCLA  360/91, 
efforts  are  being  directed  at  developing  the  necessary  graphical  routines 
and  packages  on  the  model  91  so  that  work  there  by  Center  and  LAR  per¬ 
sonnel  can  be  retransmitted  back  to  the  Center's  graphical  facilities 
attached  to  the  PDP-11,  i.e.,  the  Computek  display  scope  and  the  Gould 
electrostatic  graphics  printer.  This  effort  will  produce  systems  necessaiy 
for  activities  in  the  University’s  Atmospheric  Research  Laboratory  under 
the  direction  of  Professor  Ogura.  In  addition  to  direct  display  from  UCLA, 
a  third  party  access  through  the  Rand  Tenex  system  will  be  performed  where¬ 
by  researchers  interactively  select  their  graphics  displays  and  formats  after 
their  problem  has  been  run  and  large  output  files  produced. 
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3.  ECONOMIC  RESEARCH  GROUP 


3.1  Introduction 

During  this  period,  the  Economic  Research  Group  has  been  en¬ 
gaged  in  an  effort  to  develop  a  general  manpower,  economic  and  educa¬ 
tional  forecasting  model  which  would  be  broadly  useful  to  a  large  group 
of  users  in  many  agencies  at  the  national,  state  and  local  level,  to 
private  businesses  engaged  in  planning,  and  to  scholars  and  students  of 
economic,  manpower  and  educational  planning. 

The  Center's  econometric  research  for  the  models  used  has  been 
supported  from  several  sources,  including  the  U.S.  Department  of  Labor, 
Manpower  Administration  and  Offices  of  Policy  Development,  the  Univer¬ 
sity  of  Illinois,  and  ARPA.  Further  development  of  the  econometric 
models  will  be  supported  by  other  agencies  and  foundations.  ARPA  support 
will  be  used  to  develop  the  user  interface,  ILLIAC  IV  computational  sys¬ 
tem,  and  for  development  of  a  large  applications  system  which  is  distri¬ 
buted  over  the  Network. 

The  STEP  I  model,  which  is  almost  complete,  forms  the  first 
version  of  this  model.  STEP  II,  which  has  begun,  is  a  more  complex  ver¬ 
sion  of  STEP  I  entailing  more  sophisticated  economic  modelling  and  making 
much  greater  demands  on  computational  systems.  STEP  I  is  .implemented  on 
conventional  sequential  computers .  STEP  II  will  require  the  computational 
capacity  of  ILLIAC  IV  because  of  the  large  computational  capacity  needed 
for  matrix  operations.  The  entire  system  by  which  remote  users  will 
access  the  model  library  data  files,  and  computer  utilities  needed  for 
economic,  manpower  and  educational  planning  is  termed  MEASURE  (Manpower 
and  Economic  Analysis  System  for  Use  in  Research  and  Evaluation) .  MEA¬ 
SURE  is  to  he  developed  in  a  form  which  is  accessible  to  economists, 
planners,  budget  experts  and  others  who  are  not  trained  in  econometrics 
and  computer  science.  Properly  implemented,  MEASURE  can  provide  very 
large  economies  to  users,  and,  ultimately,  to  the  various  governments 
whose  needs  create  the  demand  for  planning  and  forecasting.  Most  impor¬ 
tant,  MEASURE  permits  the  planner  a  much  higher  degree  of  freedom  in 
developing  his  plans  since  a  much  larger  proportion  of  the  analyst’s 

4 

time  is  devoted  to  specifying  and  evaluating  his  projections  rather  than 
in  the  mechanics  of  data  collection,  model  specification,  and  data  pro¬ 
cessing. 
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It  is  by  providing  the  econometric  skeleton,  highly  user- 
oriented  computer  interface,  extensive  data  files,  and  efficient  utili¬ 
zation  of  Network  computers  ranging  from  the  PDP-10  through  the  IBM  360/91 
to  ILLIAC  IV  and  the  low-cost  communications  capability  of  the  ARPA  Net¬ 
work  that  MEASURE  will  substantially  reduce  the  costs  of  any  particular 
forecast  and  enlarge  the  community  of  planners  and  forecasters  who  can 
rise  sophisticated  forecasting  models . 

Although  ARPA  sponsorship  will  be  limited  to  system  development 
work  related  to  ILLIAC  IV,  distributed  system  experimentation  on  the  Net¬ 
work  and  user-interface  design,  we  describe  STEP  I  and  step  jj  in  some 
detail  to  provide  „he  proper  context  for  these  developments. 

3.2  STEP  I 

The  STEP  I  model  consists  of  sub-models  of  Population,  Labor 
Force,  Enrollment,  Industrial  Demand,  Demand  for  Employment  by  Industry, 
Occupational  Demand,  and  Occupational  Supply.  The  model  will  be  imple¬ 
mented  for  the  United  States  as  a -whole  and  for  the  State  of  Illinois. 

The  sub-models  are  related  as  follows:  the  Population  model 
provides  population  projections  by  single  years  of  age,  color  and  sex 
for  ary  future  year.  Hie  Labor  Force  model  uses  the  population 
projections  to  produce  labor  force  projections.  The  School  Enrollment 
model  uses  the  population  projections  to  produce  school  enrollment  by 
grade  for  any  year  in  the  future.  The  Occupational  Supply  model  uses 
the  results  of  the  Population,  Labor  Force  and  Enrollment  models  and  the 
Occupational  Mobility  model  to  generate  occupational  supply  projections 
for  any  future  year.  The  Alternative  Budget  model  generates  alternative 
patterns  of  final  demand  by  industry  for  alternative  assumptions  of  the 
pattern  of  economic  and  social  development  of  the  United  States.  The 
Industrial  Demand  model  generates  the  pattern  ;-f  demand  for  output  of 
the  various  industries  for  the  alternative  expenditure  assumptions .  The 
Employment  Demand  model  generates  the  patterns  of  employment  by  industry 
associated  with  the  alternative  industry  demand  patterns .  The  Occupa¬ 
tional  Demand  model  generates  the  alternative  patterns  of  occupational 
demand  associated  with  the  alternative  expenditure  assumptions . 

The  state  model  is  based  on  the  fundamental  assumption  that  a 
state’s  share  of  national  industry  employment  can  be  accurately  projected. 

STEP  I,  while  a  highly  simplified  version  of  reality,  repre- 


-20- 


sents  a  considerable  extension  of  forecasting  capability  for  the  user. 

It  is  the  most  comprehensive  and  general  model  of  the  employment  process 
which  has  yet  been  developed.  There  is  little  question  that  use  of  STEP 
I  in  the  preparation  of  state,  local  and  national  projections  would  re¬ 
present  a  considerable  improvement  over  existing  methods. 

The  modular  design  of  the  model  and  its  high  tolerance  for 
modified  assumptions  permit  STEP  I  to  be  easily  adapted  to  a  large  num¬ 
ber  of  special  purposes .  It  is  structured  in  a  fashion  that  permits  the 
user  a  large  degree  of  latitide  in  introducing  his  own  assumptions,  in 
extending  the  model  to  encompass  the  desired  detail,  and  to  add  on 
special  elements.  At  the  same  time,  the  model  is  itself  complete,  and 
can  be  used  for  several,  specific  applications  with  no  modification  or 
user  input  specifications  (except  to  identify  his  particular  output 
choice) . 

The  initial  demand  forecasts  prepared  using  the  demand  model 
for  the  United  States  and  for  Illinois  and  the  description  of  the  use  of 
the  STEP  I  model  will  be  published  as  CAC  technical  reports  in  the  next 
quarter. 

One  of  the  critical  problems  in  the  development  of  reliable 
forecasting  models  is  the  sensitivity  of  the  model  to  errors  of  specifi¬ 
cation.  Very  little  work  has  been  done  on  the  sensitivity  of  input- 
output  models  to  mi sspecifi cation  because  of  the  large  costs  of  invert¬ 
ing  dense  matrices  of  large  order.  The  availability  of  ILLIAC  IV  and  its 
accessibility  over  the  network  will  permit  significant  experimentation 
on  the  sensitivity  of  estimates  to  misspecifi cation  of  various  components 

of  the  model.  During  1973  the  Economic  Research  Group  will  carry  out 
/ 

extensive  sensitivity  tests  of  STEP  I  and  other  large  input-output  ma¬ 
trices  in  order  to  identify  components  of  the  model  which  are  most 
critical  to  model  forecasts. 

3.3  STEP  II 

STEP  II  represents  a  considerable  extension  of  STEP  I  in  the 
direction  of  increasing  its  realism,  and  presumably  its  accuracy  in  fore¬ 
casting  and  its  utility  as  a  research  tool.  STEP  I  relies  on  the  pro¬ 
jection  of  critical  coefficients  by  the  judgement  of  the  system  develop¬ 
ers.  The  projected  activity-industry  cr efficients,  input-output  co¬ 
efficients,  employment  coefficients,  occupational  coefficients,  birth 


rates,  mortality  rates,  labor  force  participation  rates,  and  occupational 
mobility  coefficients  are  based  on  the  judgement  in  the  development  of  the 
system.  STEP  II  envisages  forecasting  essential  parameters  of  the  model 
as  the  result  of  a  formal,  continuing  and  organized  system  of  technolo¬ 
gical  forecasting.  It  is  not  contemplated  that  the  technological  fore¬ 
casting  necessary  to  develop  best  estimates  of  the  parameters  vill  be 
conducted  in-house  at  the  Center,  but  that  the  MEASURE  system  will  pro¬ 
vide  a  vehicle  by  which  comprehensive,  consistent  and  complete  technolo¬ 
gical  forecasting  might  be  carried  out  using  STEP  II  as  a  basic  model  or 
context  for  the  forecasting .  Because  STEP  II  main!  ains  the  basic  extend¬ 
ed  modular  design  of  STEP  I,  it  trill  be  possible  for  the  technological 
forecaster  to  develop  conditional  forecasts  for  the  particular  parameters 
he  is  forecasting  without  requiring  all  other  forecasters  to  meet  and 
agree  on  basic  underlying  assumptions.  For  instance,  the  demographic 
forecaster  who  wishes  to  build  economic  influences  into  his  population 
parameters  can  specify  the  economic  preconditions  for  various  sets  of 
his  parameter  estimates. 

The  Center  plans,  however,  to  identify  teehno.logical  fore¬ 
casting  efforts  now  underway  and  to  introduce  these  groups  to  the  MEA¬ 
SURE  system.  It  is  expected  that  the  availability,  comprehensiveness 
and  econony  of  MEASURE  will  appeal  to  the  forecasters  and  induce  them  to 
join  the  MEASURE  users  group. 

In  addition  to  technological  forecasting,  STEP  II  vill  also 
differ  from  STEP  I  in  having  a  system  of  dynamic  adjustment  to  labor 
shortages  and  surpluses.  STEP  I  simply  generates  forecasts  of  occupa¬ 
tional  srq>ply  and  demand,  but  STEP  II  will  include  processes  to  predict 
the  labor  market  adjustments  which  result  from  excess  supply  or  excess 
demand. 

Moreover,  STEP  II  will  integrate. the  input-output  models  of 
flow  in  the  economy  in  the  forecast  year  with  an  annual  econometric 
forecasting  model.  STEP  I  is  based  on  uhe  distribution  of  an  aggregate 
economic  output  arrived  at  through  projection  of  aggregate  economic 
growth,  but  an  annual  econometric  model  will  restrict  the  available  eco¬ 
nometric  output  to  particular  macro-economic  aggregates,  such  as  exports, 
consumption  and  investment.  Specification  of  alternative  patterns  of 
expenditure  for  STEP  I  simply  requires  the  specification  of  alternative 
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levels,  but  STEP  II  will  require  the  specification  of  alternative  develop¬ 
mental  paths.  In  STEP  I  the  question  is, "What  will  happen  if  a  speci¬ 
fied  pattern  of  expenditures  occurs  in  1980?"  In  STEP  IJ  the  question 
is,  "What  will  happen  if  this  course  of  economic  development  occurs  be¬ 
tween  now  and  19 80?" 

The  complexity  of  STEP  II  and  its  requirements  for  large-scale 
interaction  between  experts  in  all  parts  of  the  country  should  provide  a 
major  experiment  in  the  facilitating  of  research  by  the  use  of  the  di¬ 
verse  computational  resources  and  communications  capacity  of  the  ARPA 
Network. 

3.^  MEASURE 

The  MEASURE  system  is  the  communications-software-hardware 
system  which  will  make  STEP  I,  STEP  II  and  other  econometric  models  and 
data  available  to  remote  users .  The  principal  components  of  MEASURE  are: 

(a)  The  MEASURE  operating  system  on  the  host  computer. 

(b)  The  ARPA  Network. 

(e)  The  computational  packages  distributed  over  several 
Network  computers  including  ILLIAC  IV. 

(d)  The  Model  Library. 

(e)  The  Data  Library. 

It  is  expected  that  a  MEASURE  user  will  access,  "the  MEASURE 
host  through  the  Network  to  develop  and  edit  his  particular  model  (which 
will  be  based  on  one  of  the  models  in  the  Model  Library  that  will  include 
STEP  I  and  STEP  II).  Using  his  own  data  aid  that  available  from  the 
Data  Library  the  user  will  specify  his  model.  When  he  completes  editing 
and  instructs  the  system  to  perform  the  indicated  calculations,  the  MEA¬ 
SURE  system  will  use  a  process  allocation  algorithm  to  assign  the  problem 
to  the  hose  computer,  or  a  larger  machine  (such  as  the  360/91  or  ILLIAC 
TV) .  Output  from  the  computation  may  he  stored  for  analysis  at  the 
MEASURE  host  computer.  The  development  of  the  process-allocation  algo¬ 
rithm  and  the  necessary  computational  packages  for  the  host,  360/91,  end 
ILLIAC  IV  are  principal  future  tasks  planned  under  ARPA  support.  It  is 
expected  that  the  process  allocation  algorithm  will  be  generally  useful 
in  optimizing  use  of  network  resources. 

Successful  functioning  of  the  MEASURE  system  obviously  depends 
on  the  availability  of  usable  models,  the  data  library,  and  a  community 
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of  users.  The  STEP  I  system  provides  a  considerable  advance  in  the  state 
of  the  forecasting  art  but  as  long  as  its  use  is  Hmited  to  the  local 
community,  occasional  visitors,  and  those  remote  users  able  to  pay  the 
substantial  costs  of  normal  long-distance  telephone  connection  it  will 
not  be  fully  utilized.  There  is  little  benefit  in  developing  the  soft¬ 
ware  for  such  complex  systems  unless  the  essential  components  are  likely 
to  become  widely  available  and  unless  there  is  effective  demand  for  the 
system.  At  the  same  time,  it  is  difficult  to  measure  effective  demand 
for  such  a  system  until  it  is  specified,  vith  a  firm  date  for  avail¬ 
ability  and  narrowly  defined  user  costs. 

In  1972-73  the  Center  will  be  engaged  in  experiments  with 
several  remote  users  of  the  STEP  I  system.  These  include  academic  re¬ 
searchers  at  other  universities  and  government  agencies.  The  experi¬ 
ence  gained  during  this  year  will  aid  us  in  gaining  understanding  of  the 
problems  of  using  large,  complex  systems  in  a  research  mode  over  the  ARPA 
Network . 
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ILLIAC  IV  LANGUAGE  DEVELOPMENT  FOR  THE  PHASE  II  SYSTEM 


During  the  latter  part  of  this  period,  a  study  began  to  address 
the  need  for  a  convenient  linguistic  and  systemic  means  to  specify  and 
control  timely  flow  of  information  between  hierarchic  memory  levels  in 
support  of  a  non-core  contained  computation  on  ILLIAC  IV.  Incorporation 
of  run  time  hardware  allocation  constraints  into  the  data  flow  process 
is  also  a  problem. 

As  an  ideal  solution,  one  would  like  to  have  a  compiler  which 
would  accept  the  source  program  and  the  run-time  storage  allocations,  and 
produce  the  appropriate  set  of  ILLIAC  IV,  Memory  Management  Processor 
(MMP),  PDP-10  and  CCA  datacomputer  "programs"  necessary  to  efficiently 
support ,  synchronize  and  perform  the  indicated  computation.  Work  during 
the  next  six  months  will  attempt  to  identify  the  framework  and  method¬ 
ology  necessary  to  support  attainment  of  these  goals. 


5.  ILLIAC  IV  INFORMATION  RETRIEVAL  AliD  STATISTICi\L  SYSTEM  (IRSS) 


5.1  Summary 

An  investigation  was  made  and  completed  on  information  re¬ 
trieval  using  ILLIAC  IV.  The  conclusions  reached  were  that  the  machine 
is  not  sui'.  ed  to  the  general  problems  of  information  retrieval  and  data 
management  due  to*  the  awkward  I/O  structure,  i.e.,  the  lack  of  parallel¬ 
ism  in  I/C.  However,  files  that  can  be  linearized  for  serial  searches 
are  applicable  to  the  ILLIAC  since  PE  parallelism  can  be  exploited.  Ser¬ 
ial  searching  is  no  longer  a  general  problem  of  data  management  due  to 
random  access  equipment  and  techniques.  If  serial  search  is  necessary, 
it  can  be  done  on  the  TV;  however,  its  use  is  probably  limited. 

In  November,  functional  and  detailed  design  of  the  IRSS  began. 
Progress  is  slow  as  the  SSystem  has  dependencies  on  the  datacomputer  and 
the  Ames  Phase  II  System  which  are  not  completely  defined. 

The  scope  of  the  IRSS  has  been  decided  and  includes  the  follow¬ 


ing: 

(a)  The  IRSS  will  relate  to  information  retrieval, 
not  to  information  management. 

(b)  It  will  provide  a  user-language  (non-procedural) 
and  recognize  data-item  names. 

(c)  The  IRSS  will  provide  data  manipulating  routines 
(statistical)  identifiable  through  the  user- 
language,  operable  in  the  ILLIAC  IV. 

(d)  The  IRSS  will  he  developed  with  the  Ames  Phase 
II  system  as  its  base. 

(e)  The  datacomputer  will  bs  used. 

(f)  The  IRSS  will  provide  a  means  for  researchers 
whose  data  reduction  problems  are  large  enough 
to  warrant  use  of  the  ILLIAC  IV  to  gain  access 
to  it. 

(g)  Data  transposition  routines  will  he  available 
for  users  who  need  a  specific  format  of  data 
in  PE  Memory. 

(h)  The  IRSS  will  he  interactive  but  will  not 
guarantee  response  in  terms  of  completion. 

-26- 


With  tliis  emphasis  ,  the  next  period  will  involve  the  develop¬ 
ment  of  the  user-interface  that  is  the  non-procedural  language.  Based  on 
the  current  state  of  the  Ames  I/O  design,  detailed  design  of  the  inter¬ 
face  to  the  ILLIAC  IV  operating  environment  will  he  deferred.  Design  and 
implementation  of  the  language,  symbol  tables,  parser,  and  translation  to 
data  language  will  be  investigated. 

An  initial  version  of  a  user-oriented  interactive,  file  hand¬ 
ling  system  called  MONICA  was  completed  during  this  reporting  period. 
MONICA  was  developed  to  provide  researchers  easy  access  to  a  set  of  statis¬ 
tical  subroutines.  Emphasis  was  given  to  the  user-interface  and  this 
experience  will  be  used  in  developing  the  IRSS  user-interface.  A  user's 
manual  for  MONICA  will  be  produced  in  the  next  quarter. 

Work  was  begun  on  defining  user  needs  for  an  ILLIAC  IV  Statis¬ 
tical  System  and  for  reviewing  existing  systems.  The  ILLIAC  Statistical 
System  will  be  an  integral  part  of  the  ILLIAC  Information  Retrieval  Sys- ' 
tern. 

Dr.  David  Beckles,  presently  of  the  University  of  Nevcastle- 
on-iyne,  England,  will  join  the  staff  of  the  Center  in  August  to  take 
primary  responsibility  for  the  design  and  implementation  of  the  ILLIAC 
Statistical  System.  In  January  of  this  year,  discussions  were  held  with 
a  number  of  eminent  users  and  designers  of  large  statistical  systems. 

The  consensus  of  this  series  of  meetings  was  that  it  would  be  unwise  to 
implement  a  national  resource  like  the  ILLIAC  System  without  continuous 
consultation  with  the  user  community.  We  intend,  therefore,  to  set  up 
from  the  user  community  a  statistical  advisory  committee  to  provide 
comments  and  suggestions  to  the  design  group.  Dr.  John  Sonquist,  Pro¬ 
fessor  of  Sociology  at  the  University  of  California  at  Santa  Barbara  has 
agreed  to  head  the  advisory  committee. 
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6.  ILLIAC  IV  IMAGE  PROCESSING  SYSTEM 


6.1  Summary 


At  the  request  of  NASA’s  Earth  Observations  Program,  the  Cen¬ 
ter  is  presently  exploring  the  application  of  the  ILLIAC  IV  hardware  sys- 
J  -i  to  the  large-scale  data  storage  and  processing  problems  associated 
the  automatic  interpretation  and  management  of  high-altitude  air- 
c.  t  and  satellite  earth  resources  imagery . 

As  background  for  this  effort,  the  Center  has  surveyed  exist¬ 
ing  methods  of  earch  resources  data  collection  by  remote  sensing,  current¬ 
ly  operational  hardware  and  software  systems  for  data  processing,  inter¬ 
pretation  and  information  management,  and  meaningful  applications  of  the 
technologies  to  specific  problems  of  social  concern.  Preliminary  inves¬ 
tigations  have  shown  the  ILLIAC  IV  capable  of  efficiently  processing,  in¬ 
terpreting  and  managing  large  quantities  of  multi-spectral,  digital  imagery. 
The  ARPA  Network  seems  a  practical  means  for  interfacing  the  numerous, 
geographically-dispersed,  potential  users  of  such  a  system. 

As  a  result  of  these  initial  determinations,  the  Center  is  pre¬ 
paring  a  proposal  to  NASA  (and  perhaps  to  other  federal  agencies)  for  the 
design  and  implementation  of  a  prototype  ILLIAC  IV-based  information  pro¬ 
cessing  system  to  support  NASA's  Earth  Observations  Program.  The  Labora¬ 
tory  for  the  Applications  of  Remote  Sensing  of  Purdue  University,  supported 
by  NASA,  is  assisting  the  Center  in  this  effort  and  has  agreed  to  lend 
methodological  expertise  to  future  efforts .  NASA-Ames  Research  Center 
has  also  participated  in  preliminary  explorations  and  has  expressed  sub¬ 
stantial  interest  in  such  a  proposal. 
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7*  EDUCATIONAL  ACTIVITIES 

7.1  Documentation 

(a)  Volume  1  of  "An  Introductory  Description  of  the  ILLIAC 
IV  System"  was  published  on  July  15,  1971.  Two  hundred  copies  were  sent 
to  Ames  and  three  hundred  copies  were  kept  at  the  University  of  Illinois. 
Following  is  an  abstract  of  the  hook: 

Written  specifically  for  an  applications  programmer,  the  hook 
presents  a  tutorial  description  of  the  ILLIAC  IV  System.  Volume  1  con¬ 
tains  three  chapters-Background,  Hardware  Structure  and  the  Assembly  Lan¬ 
guage,  ASK,  as  well  as  a  Hardware  Glossary.  Mary  illustrative  problems 
are  used  to  educate  the  beginner  in  the  use  of  the  ILLIAC  IV  System. 

(b)  A  very  short  "Ten  Page  Description  of  the  ILLIAC  IV  Sys¬ 
tem"  is  now  in  preparation.  It  was  synopsized  from  "ILLIAC  IV  Hardware 
Structure,"  presented  at  the  Ninth  Annual  Allerton  Conference  on  Circuit 
and  System  Theory  in  October  1971  and  from  "The  ILLIAC  IV  System," 
written  for  the  IEEE. 

7-2  Teaching  and  Consulting 

(a)  CS  1*91,  "Architecture,  Applications  and  Languages  for  a 
Parallel  Computer"  was  offered  at  the  University  of  Illinois  as  a  graduate 
course  in  the  Computer  Science  Department.  The  course  has  now  been 
changed  to  CS  397  "The  ILLIAC  IV  Computer  System"  and  is  being  taught 
this  semester. 

(b)  Short  ( one-hour  to  one-day),  informal  seminars  dealing 
with  the  ILLIAC  TV  System  were  given  to  University  of  Illinois  staff  and 
students.  Some  consulting  concerning  ILLIAC  IV  was  done  by  telephone 
with  people  outside  the  University  community. 

7.3  Training 

Two  Affirmative  Action  programmer  trainees  are  being  trained 
in  the  arts  and  crafts  of  programming  for  a  digital  computer  and  in  algo¬ 
rithmic  problem  solving. 
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